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O 

We analyse the bound states of negatively charged ions which were predicted to exist 
because of the presence of a magnetic field by Avron et al fjj]. We confirm that the 
number of such states is infinite in the approximation of an infinitely heavy nucleus and 
, provide insight into the underlying physical picture by means of a combined adiabatic and 

(■ — | perturbation theoretical approach. We also calculate the corresponding binding energies 

which are qualitatively different for the states with vanishing and non-vanishing angular 
■ momentum. An outlook on the case of including center of mass effects is presented. 

On: 



^ 1 I. INTRODUCTION 

o 

The behaviour and properties of negative ions became during the past years a branch of intense research. 
. There exists an enormous interest in the electronic structure and dynamics of negative ions, both from 
the theoretical as well as experimental point of view (see, e.g., the review []2| and the references therein). 
According to our present knowledge about atomic ions it is most likely that singly charged negative 
ions possess in the absence of a magnetic field only one stable ground-state configuration. For example 

^sO ' for the H~ ion, this state is the l 1 ^ electronic state, and a rigorous proof that this the only possible 
bound state was given in ||]. Also, some atoms, like for example Be, N, Ne, Mg, Ar do not possess any 
stable negative ion state (see, e.g., ^|||). On the other hand, one can expect that in the presence of a 
magnetic field a lot of new discrete energy states of negative ions can appear. This expectation is based 
on the statement that for any negatively charged ion the number of discrete energy states is infinite in 
the presence of a magnetic field [0. However, this statement was formulated as the conclusion of a formal 

• i-h ' mathematical treatment which does not provide a transparent physical picture of the appearance of the 
infinite sequence of bound states nor does it estimate the corresponding energies. Although a wide set 

^; 

Oh 



X 



of works was focused on the influence of the magnetic field on low- lying ion states (e.g., J6|-|l2|), there is 
no systematic treatment of the highly excited anions predicted in |jj . In the present paper we develop a 
physical approach showing transparently from which quantum mechanical grounds the infinite sequence 
of bound anion states appears in the presence of a magnetic field. In Section 2, we first reduce the 
problem to a one-particle, and then to a one-dimensional one, where to analyse the bonding one has to 
consider the motion of the external electron along the magnetic field in an effective potential depending 
on the quantum state. We confirm that, in the approximation of the static (infinitely heavy) nucleus, 
the number of the bound states is infinite, and derive estimates for the corresponding binding energies 
valid for arbitrary ions. In Section 3, we apply these estimates to the ion H~ and verify them treating 
the one-dimensional motion of the external electron numerically. We also discuss the bound states of 
an exotic ion formed by attaching the muon to the hydrogen atom. Besides, we analyse the importance 
of the non-adiabatic coupling between the states of the external and atomic (core) electrons. A general 
discussion and concluding remarks are given in Section 4. 
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II. GENERAL 



It is well-known (see, e.g., jl 13 1) that the electron affinities for the negative ions are usually much 
smaller than the binding energies of the corresponding neutrals. This fact implies that the extra electron 
is weakly bound to the atom and its behaviour and properties strongly differ from that of internal (core) 
electrons. For the H~— ion, the binding of the extra electron is exclusively due to correlation between 
the two electrons. In the following we will study the highly excited states of the negative ions which 
appear in the presence of a magnetic field and which, as we will estimate, correspond to very small 
electron affinities. We therefore use the approximation of a weakly bound external electron neglecting its 
exchange interaction with the core electrons. To exploit such a model, it is convenient to split the total 
Hamiltonian for a singly negatively charged ion into three terms, 



Here the first term, 



H ion = H at +H C + W. (1) 
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a=l 



describes a neutral atom with infinitely heavy nucleus in a magnetic field B with the vector potential 
A a = A(r a ). Summation is carried out over all the atomic (core) electrons labelled with the subscript 
"a", Z is the nucleus charge number, m G and — e = — |e| < are the electron mass and charge. The 
potential energy U includes the Coulomb energies of the interaction of the core electrons with each other 
and with the nucleus. S is the total spin of the atom, and /xb = eh/(2m e c) is the Bohr magneton. The 
second term in (Q), 

1 



2m r 



c 



p+-A +2^ B Ber, (3) 



corresponds to an extra (with respect to the atom) electron, er is the electron spin. The last term, 

^ r - r a r 

describes the Coulomb coupling of the extra electron to the atom. 

When the extra electron is weakly bound to the atom, the character of its motion differs strongly from 
that of the core electrons. In particular, the external electron can be assumed to move much slower 
than the core ones. In this case, the ion states can be successfully described in terms of a quasimolecular 
approach, and the problem of binding can even be reduced to a one-particle one. Below we briefly describe 
this approach. 



A. Quasimolecular approach to the problem of binding 

Let us consider the Hamiltonian 

H' at = H at + W. (5) 

It does not include the kinetic energy of the external electron, therefore the latter can be considered as 
static in space, so H' at describes the motion of the core electrons only. Since the potential energy term 
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of H' at depends parametrically on the position r of the external electron, the energy spectrum of this 
Hamiltonian as well as the corresponding eigenfunctions also depend on r. Let us denote them by 

EKr) = (i\H'Ji) (6) 

and 0»(ri, r 2 , . . . , r z ; r), respectively, where the index i labels the eigenstates of the Hamiltonian H' at 
with some choice of r. Since at any r the states \i) compose a complete basis set, an eigenfunction of the 
total Hamiltonian (Q) can be presented as 

*(n, r 2 , . . . ,r z ;r) = Y^ A' (r)(t>i> (n, r 2 , . . . , r z ; r), (7) 

with the expansion coefficients i^k' ( r ) depending on the position of the external electron. It should be 
noted that this wave function is not completely antisymmetric with respect to all the ion electrons, only 
its "atomic part" , 4>i' can be assumed to be properly antisymmetric with respect to the core electrons. 
Complete antisymmetrisation would significantly complicate the further consideration. On the other 
hand, the main yield of such a complication would be taking into account exchange interaction between 
the external and core electrons which is definitely negligible for states for which the external electron is 
localized far from the atom. Notice also that the form (|7|) of the ion wave function is exact if another 
charged particle different from the electron (like for example a muon) is attached to the atom. 

Substituting the wave function (ffl) into the Schrodinger equation with the Hamiltonian (0), subse- 
quently multiplying from the left-hand side by (f>* and integrating over the positions of the core electrons 
we will arrive at the set of coupled equations 

[H e + E[{r) + h u (r) - E tot ] <fc(r) = - ^ hu> </v(r). (8) 

Here l?tot is the total eigenenergy of the system, and the matrix elements hi i> are given by the following 
equation 

1 -<&V<M + 7^-v#Ni'} + — (<Pi\P<Pi')P , (9) 



2m c Zm c c m c 

where I = r x p is the angular momentum of the external electron. Notice that if both i and i' relate to 
the ground state of the atom, |0), then the last two terms in Eq. (|J) equal zero while the first term can 
be transformed as follows 

2 



R ' (10) 



^(j}{ri,r 2 , ...,r z ;r) 



In the following we restrict ourself by considering the attachment of the extra electron to the neutral 
atom in the ground state. In this case, the right-hand side of Eq. (||) is associated with the coupling 
between the ground and excited states of the core electrons due to the motion of the external electron. 
Assuming the ground state of the neutral atom to be energetically well-separated from the excited states 
this coupling can be neglected. Then Eq. (|[) becomes essentially a one-particle Schrodinger equation, 

[H e + V (r) + hoo(r) - E ] Mr) = 0, (11) 

where, for the sake of convenience, we have shifted the energies by the ground state energy of the isolated 
atom, E^ = E' (r — » oo), and introduced 

E Q = - 4 0) (12) 



3 



and 



V (r) = E' (r) - E^. (13) 

Equation (^lj) describes the motion of the external electron in the magnetic field and the potential which 
consists of two parts: the static term, Vo(r), and the non-adiabatic (dynamic) correction, hoo(r). The 
latter, as we shall directly estimate for the H~— ion, decreases with increasing r faster than V{r) and can 
be neglected at distances exceeding the atomic size. Therefore, to analyse the binding mechanism for the 
excited ion states one has to solve the Schrodinger equation 

[H e + V(r) - E] Hr) = 0, (14) 

where for brevity we have omitted the index "0" related to the ground state of the core electrons. 

At distances of the extra electron from the atom strongly exceeding the atom size, one can consider 
the operator (Q) as a perturbation of the atom and evaluate it by the multipole expansion, 

* — + cn> 

where n = r/r is the unit vector in the direction of r and the indices a and (3 run over the Cartesian 
coordinates. 



D = -eJ2r a (16) 

o=l 

and 

z 

Qa/3 = -ej^ (rlSa/3 - 3x aa x a p) (17) 

a=l 

are the operators of the dipole and quadrupole momenta of the atom, respectively, x aa and x a p denote 
the components of r a , r a — (x a i, x a 2, %a,3)- Then the potential in Eq. ( fl4"|) can be approximated as 

v , \ _ e(D)n e{Q af3 )n a n j3 e 2 K a pn a np 

where (D) and (Q a p) are the mean values of the dipole and quadrupole momenta, respectively, of the 
unperturbed atom described by the Hamiltonian (|^), and n a p is the polarizability of the atom in an 
electric field. The first two terms in ( |l8| ) represent the first-order perturbation corrections. The last term 
in (fL8|) is the second-order correction, it corresponds to the dipole term in (^5j) treated as the perturbation 
of the atom by the electric field £ = —en/r 2 . 



Expression ( jlq ) can be further simplified if we take into account that the potential (13) has its symmetry 
axes directed along the magnetic field. A simple argumentation of this statement is the following. In 
the presence of the magnetic field it is natural to assume that the axis of quantization for the atom is 
directed along B. Then the potential (|l^) depends parametrically on the two vectors, B and r. Since 
the potential is scalar it must depend only on scalars which can be constructed from these vectors. These 
scalars are B 2 , r 2 and Br. This means that if B is directed along the z— axes, then V(r) — V(r±,z). 
As a result the non-zero components of the mean dipole and quadrupole momenta of the neutral atom 
can only be (D z ), (Q xx ) — (Qyy) an d (Qzz), respectively. Also, the polarizability takes the diagonal 
form with the components n xx = K yy = k± and n zz — Kii. Furthermore, since the z— parity of the core 
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electrons is the integral of motion when they are not perturbed by the interaction with the external 
electron, we have (D z ) — 0. Then Eq. (|l^) can be transformed as follows: 

e 2 X(2z 2 -r 2 ) e 2 
V(r ±7 z) = V ^ ±J - (k± sin 2 + K „ cos 2 0) , (19) 

where ■& is the angle between n and B and the coefficient A can be expressed in terms of the mean values 
of the squares of the longitudinal and transverse coordinates of the core electrons of the unperturbed 
atom, 

X = j^[2{zl)-{rl x )] . (20) 

a=l 



B. Binding in a magnetic field as ID problem 

The Hamiltonian that determines the one-particle Schrodinger equation ( |l4| ) explicitly reads 



H eS = + ^L. + V(r ± ,z), (21) 



where 



7v ± =p ± + ^-Bxr ± (22) 
Zc 

is the transverse kinetic momentum of the electron. We have introduced the symmetric gauge of the 
vector potential, A = (l/2)Bxr, being the most appropriate one because of the axial symmetry of the 
potential V. We also have omitted the spin part of the Hamiltonian which determines the trivial shift of 
the ion energy spectrum and does not affect the binding energies. 

Since the potential in the Hamiltonian (|l]) decreases rapidly with increasing distance of the extra 
electron from the atom, even quite a weak magnetic field can influence the transverse motion of the 
extra electron to a much larger extent than the atomic potential. In this case the eigenfunction of the 
Hamiltonian (|2l| ) can be efficiently expanded in terms of the the Landau states, (T"j_|n, s), of the motion 
of the electron across the magnetic field. These states are the common eigenstates of the operators tt]_ 
and l z , where I = rxp is the electron angular momentum, and labeled with the two quantum numbers, 
n and s, that determine the corresponding eigenvalues, 

tt^ = (?i/a m ) 2 (2n + l) , n = 0,l,2,..., 

l z = —hs , s — — n, —n + 1, — n + 2, . . . , (23) 



a m = yj cfi/(eB) is the magnetic length. The first eigenvalue determines the Landau energy spectrum 

^ n = ^ = ^M) • (24) 



2m c m c c \ 2 / 

Hence n is called the electron Landau level number. 

Because of the axial symmetry of the atomic potential the longitudinal component of the electron 
angular momentum is an integral of motion for the Hamiltonian (pi]), [l Zl H c s\ = 0, and the quantum 
number s can be used to label the eigenstates of H c g, i.e. E — E s . The expansion of the corresponding 
eigenfunctions over the Landau states thus involves only different Landau level numbers and reads 
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^Ps(r) = ^2g ns {z){r ± \n,s). (25) 

a 

The expansion coefficients g n s{z) and the energy E can be found by solving the system of coupled 
equations, 



,,2 



^- + V^(z) + E^-E s 



g ns (z) = -J2Vn S rhz)9n> s (z), 11 = 0,1,2,..., (26) 



where 



V$(z) = (n,s\V(r ± ,z)\n',s) (27) 



are the effective longitudinal potentials obtained as the matrix elements of the atomic potential V(r±, z) 
within the space of the Landau states of the external electron. 

The quantity i?Q an = Tjws/2, where uib = eB/(m c c) is the electron cyclotron frequency, determines the 
continuum threshold for the states of the external electron. The states of this electron belong to different 
manifolds associated with the different Landau energies £^ an , n = 0, 1, 2, . . .. The bound states can be 
classified by the numbers (s, n, v) where v equals the number of nodes of the function g n s(z) for the 
leading term in the expansion ( p5| ) - such a classification is similar to that used in treating hydrogen-like 
atoms in a strong magnetic field, see, e.g., jT3. 

We can expect (and our results for the ion H~ confirm this), that the binding energies for the ionic 
states induced by the presence of the magnetic field are small compared to the electron cyclotron energy, 
Hlub- Therefore, all such states associated with the manifolds n = 1,2, ... lie in the continuum, i.e. are 
not bound. We thus will focus on the states related to the ground Landau manifold, n — 0. While 
treating them we can further neglect the coupling of the ground Landau manifold to the higher ones, 
e.g., omit the sum on the right-hand side of Eq. (Eq). This approach yields the ID Schrodinger equation, 



g s (z) = 0, (28) 



which describes the motion of the external electron along the magnetic field. In this equation, we have 
shifted the energy by the zero-point Landau energy, 

e s =E S - E^ n . (29) 

Then, when e s is negative, the value of — e s is the binding energy of the external electron. Also, for brevity 
we have omitted the index n = and define the effective longitudinal potential by V s (z) — Vqq\z). In 
order to get an idea of the binding properties of this potential one can use the weak coupling one- 
dimensional theory (see, e.g., |lj]). This theory says that if the potential vanishes at \z\ — > oo and the 
integral 



I s = / V s (z)dz (30) 

J — oo 

is negative than there is at least one bound state of finite longitudinal motion, and the estimate of the 
corresponding binding energy for that state is 

-s s =mJ 2 J{2h 2 ) . (31) 

As we shall see, the condition I s < holds for infinitely many possible quantum numbers s and thus the 
number of bound states associated with different values of s is infinite, in accordance with the conclusion 



G 



of M. For each s, the estimate (|3l]) gives the binding energy for the state for which the longitudinal part 
of the wave function of the external electron has no nodes, e.g. v = 0. 

To exploit the latter analytical approach let us introduce the transverse probability density for the 
states related to the ground Landau manifold 



^=jy^ n =°^=MM^(-tk) ■ (32) 

where the integration is done over the electron azimuthal angle. In terms of the density ( |32"| ) and the 
atomic potential V(r±,z) we have 
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I s = r ± dr ±Ps (r ± ) dzV(r ± ,z). (33) 



o 

Let us now adopt an approximation of the atomic potential given by Eq. ( |l9| ) for large distances. It is 
easy to check that at any r± / 

(2z 2 - r 2 , ) dz 

HjT- = 34 

(rl + z 2 f 2 



and thus the first, quadrupole, term in (19) does not contribute into I s . This implies that the polarization 
part of the atomic potential plays a key role in binding the external electron. Let us further assume that 
the influence of the magnetic field on the atomic polarizability is negligible. As it is known fig] , in 
the absence of the magnetic field an atom in the ground state of the zero total angular momentum has 
an isotropic polarizability, k, i.e. we assume k± ps «|| « k. Then the polarization term in ( |l9| ) is net 
attractive and the condition I s < is fulfilled. Straightforward integration in Eq. (J33|) and use of the 



relation (31) gives 

-z* = ww r VV Ryl ' (35) 

where T is the complete gamma-function, Ry = m c e 4 /(2h 2 ), 7 = a^/a^ = B/Bo, and Bq = 2.35 x 10 5 T. 
This estimate is not appropriate for s — because of the divergence of the integrations at the origin 
r± = z = for the potential ([II]). However, for s > above value of the binding energy is a good 
approximation for 7«1 (laboratory magnetic fields) since the main contribution to the integral comes 
from distances sufficiently larger than ao, the natural measure of the size of the neutral atom. Using 
atomic units (2Ry for the energy and Bq for the magnetic field strength) we rewrite Eq. ( |35| ) in a recurrent 
way, 

- e s = 0.1211k 2 B 3 S 2 , s=l,2,..., (36) 

where 

Si = 1 , 

S s = [1 - (1.5/*)]*._i , s = 2,3,... . (37) 

One should notice that the scaling behaviour of the binding energies with the magnetic field strength, 
—e s cx -B 3 , coincides with that predicted in M. Also, at large s we have S s « s~ 3 / 2 and thus obtain the 
behaviour of the binding energies for large s, — e s oc s~ 3 . Such a sharp decrease of the binding energy with 
the quantum number reflects the fact that in the plane perpendicular to the magnetic field the external 
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electron follows the Landau orbit with the expectation value for the radius squared (r^) = 2(s + l)a^, 
and at large s the transverse probability density ( p2| ) becomes a peaked function of r\ centered around 
r\_ = (r 2 ^). Therefore, as s increases the external electron is bound to the atom at progressively larger 
distances from it. 

We can also establish the behaviour of the binding energy with the magnetic field strength for the 
state s = 0. In this case, contrary to the limit of large s, the transverse probability density varies very 
slowly with r± and the scale on which it reduces significantly from its maximal value at r±_ — is the 
magnetic length, a m , which strongly exceeds the atomic potential scale, clq. Therefore, when calculating 
the integral ( |33| ) we can replace the transverse probability density by its maximal value, p s (0) = 
In this way we obtain 



r /■ \ / 'x ~ 



£0= 2 



+00 



dr± / AzV(r^,z) 



B 2 , (38) 



where we have used atomic units again as in the remaining part of our work. We remark that the domain 
of the neutral atomic core contributes quite significantly to the integrations in ( |38| ) and thus the use 
of the quasimolecular (adiabatic) picture to describe the motion of the extra electron becomes a crude 
approximation. However, a more accurate treatment is expected to change only the numerical coefficient 
in Eq. (^) . The scaling behaviour of the binding energy with respect to the field strength is expected to 
be the same, i.e. — £0 oc B 2 . 



III. BINDING ENERGIES OF THE H AND H/x" IONS 

The approach developed in the previous Section can be directly applied to the H ion. With a good 
accuracy we can neglect the influence of the magnetic field on the ground state of the core electron, and 
use the well-known value (see, e.g., JlJ]), n = 9/2, for the polarizability of such a configuration. Then 
the estimate ( |36| ) gives 

-e s = 2A53B 3 6 2 . (39) 

More reliable results for the binding energies for the states of the electron attached to the hydrogen 
atom which are induced by the presence of the magnetic field, including the state s = 0, can be obtained 
if one knows the atomic potential V. When the influence of the magnetic field on the hydrogen atom is 
neglected this potential is spherically symmetric, V — V(r), and can be found in the literature. We use 
the results of Walles, Herman and Milnes |l^] who rigorously calculated the energy levels of an electron 
moving in the field of two fixed, at distance r from each other, charges +e and — e using the separability 
of the problem in confocal elliptic coordinates. Their electronic ground state energy, -EwhMi is directly 
related to our potential V(r), 

K(r)=£ W HM(r)-- + 0.5, (40) 
r 

when adding the Coulomb interaction energy between the two static charges and the hydrogenic ground- 
state binding energy (cf. Eq. (fllf)). The energy -Ewhm was tabulated in the wide range of r, 0.840380 < 
r < 30.0. In Figure 1, we show the corresponding values of V by dots and use the spline interpolation to 
determine the value of the potential between them (solid line) . At larger r the potential can be evaluated 
by its polarization tail, 
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while at the lower r one can put Ewhm = in Eq. (|4^) and get 

V(r) = --+0.5 . (42) 
r 

As Figure 1 demonstrates, these extrapolations match very good the data from Jl8| . 

With the well-defined potential V(r) we can find the value of the coefficient in the estimate of the 



binding energy for the state s = 0. Numerical integration of the potential in Eq. (38), with the natural 
replacement 2irr±dr±dz — > 4irr 2 dr, yields 

- s = 6MB 2 . (43) 



To verify the estimates (39) and (|4^) we have also solved numerically the Schrodinger equation (pq), 
calculating the effective potentials V s (z) numerically from the potential V(r) presented in Figure 1. First 
of all, we were searching for the states with no nodes of the longitudinal wave function of the external 
electron (y — 0— states). The results for the corresponding binding energies, for different numbers of s 
and field strengths, are shown in Figure 2 by open circles. Solid lines represent the estimates ( |39"| ) and 
([43|), and we can conclude that they are in fairly good agreement with the numerical results. We have also 
performed the search for the bound states with higher longitudinal excitations, v = 1,2,..., but found no 
one, for the magnetic field strengths 1CP 6 < B < 10~ 4 . This indicates that the effective potential V s (z) 
is so weak that for each s it can bind the external electron only in the v = 0— state. 

Finally, to control the validity of the quasimolecular approach, we have calculated the quantity h(r) = 
hoo{r), the non-adiabatic correction to the potential V(r), determined by Eq. (|Io|). For this purpose, the 
wave functions of the electron moving in the field of two fixed charges were computed via a numerical 
two-dimensional grid method which was invented by M.V. Ivanov (for details of the method as well as 
the description of the corresponding code ATMOLMESH we refer the reader to p9| ). The corresponding 
results are shown in Figure 1 by open circles. To confirm that the numerical results reproduce correctly 
the scaling behaviour of h with r, we also give the perturbation estimate, 

h(r) = A , (44) 

whose derivation is outlined in the Appendix. With increasing distances exceeding the size of the neutral 
hydrogen atom, the non-adiabatic correction becomes smaller and smaller compared to \ V (r) \ . One should 
notice that at small r the non-adiabatic correction exceeds the value of |V(r)|. This implies that the 
non-adiabatic coupling effects are quite significant for the state s — 0. However, since the corresponding 
binding energy obtained neglecting these effects is very small, we believe that a more accurate approach 
would only change the numerical coefficient in Eq. ( ff3|) but not the scaling low, — So oc B 2 . 

The approach developed in this paper can also be directly applied to the analysis of highly excited 
exotic anions, for example to the ion formed by attaching the muon to the hydrogen atom. Due to the 
fact that the shell of the Hfi~— ion is formed by two different particles (the electron and muon) there is 
no need for antisimmetrization (exchange) and the ansatz (^) for the wave function of the ion becomes 
exact. Moreover, since the muon is much heavier than the electron (the muon-to-electron mass ratio is 
"V/m c ~ 207), the binding energies for the muonic ion are larger. Replacement of the electron mass by 
the muon mass in Eq. ([}l]) leads to the following estimates 
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- £ = 1306B 2 , s = , 
-e s = 508B 3 S 2 s , a = 1,2, 



(45) 



In Figure 2, the corresponding energies are plotted as dashed lines. Evidently, the non-adiabatic correction 
to the one-particle potential for the Hyx - ion is m^/m c times smaller than that for the ion H~ (see 
Figure 1), which makes the estimate for s = more reliable in the case of the H/i~ ion. 

IV. SUMMARY AND OUTLOOK 

In this paper, we have studied the states of the negative atomic ions which are induced by the presence 
of a magnetic field. Apart from the neglect of the exchange interaction between the external and atomic 
(core) electrons the equations we start with are exact. Performing an analysis of the weakly bound states 
of the external electron, it is possible to reduce the problem to an effective one-particle by neglecting 
the non-adiabatic coupling terms. In the presence of the magnetic field which determines the transverse 
motion of the external electron further simplifications are possible which reduce the question of the 
binding mechanism to a one-dimensional Schrodinger equation for the motion of the external electron 
along the field. 

We can conclude that in the approximation of the infinitely heavy nucleus the number of the ionic bound 
states is infinite, in accordance with the general theorem of [Q. A quite simple and appealing physical 
picture of the appearance of such states is that the external electron can be attached to the atom with 
different values of its angular momentum along the magnetic field. The number of such possible definite 
values is infinite and for each value a different ID effective potential appears. This potential determines 
the motion of the external electron along the field and can bind the electron in at least one quantum state. 
We have labelled these states by the integer number s = 0,1,2,... (negative of the magnetic quantum 
number of the external electron), and obtained general estimates of the binding energies of negatively 
charged ions. These estimates establish both the dependence of the binding energies on the magnetic 
field strength and on the quantum number s: for s = the binding energy scales with B as B 2 and for 
s = 1, 2, . . . the scaling low is — e s oc B 3 ; at large s the binding energies behave as — e s oc s~ 3 . To apply 
the estimates for s > to a specific ion one needs only to know the polarizability of the corresponding 
neutral atom in an external electric field, while the estimate for s — involves the whole atomic potential 
acting on a static external electron. 

We have applied the estimates obtained to the H~— ion, exploiting the known polarizability of the 
hydrogen atom and the potential for the static two-charge Coulomb problem (jl8). We also have verified 
that these estimates are in complete agreement with results of the numerical integration of the Schrodinger 
equation for the motion of the external electron along the magnetic field. 

The binding energies for the H~— ion that we have obtained are indeed very small compared to the 
binding energy of the hydrogen atom. At the largest magnetic fields available at laboratories now, 
B ~ 30 T, the binding energy for the H~— ion is approximately 2.8 x 10~ 3 mcV for the state with the 
zero angular momentum of the external electron, and it is « 1.4 x 10 -7 meV for the state with s = 1. 
The magnetically induced bound states of the exotic muon ions possess a much larger binding energy. 
For the same value of the magnetic field strength, the muon affinity to the hydrogen atom is 0.58 meV 
for the state s = and 2.9 x 10~ 5 meV for the state s = 1. 

It also looks quite challenging to detect the excited anion states for heavier atoms which have polar- 
izabilities significantly larger than the hydrogen atom. For example, the polarizability of Cs is, with 2% 
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accuracy, k = 403 |T^|. The resulting electron affinity to this atom for the state s = 1 at B — 30 T is 
k, 2 x 10~ 4 meV. Another challenging example is the anion of Ba, which does not exist at B = 0. The 
polarizability of Ba, with 8% accuracy, is k = 268 J!?]] , and the binding energy of its excited anion, for 
the external electron state with s = 1 at B = 30 T is ~ 1.8 x 10 -5 meV. The corresponding states of 
the exotic ions, Cs/i~ and Ba,/j,~ , possess the binding energies of « 4 x 10~ 2 meV and «4x 10~ 3 meV, 
respectively. 

Although we cannot directly apply the estimate (|38| ) for the binding energy of the external electron 
with zero angular momentum to heavier ions, in order to have an idea about their order of magnitude we 
can speculate that these binding energies differ from the result ( [43| ) by the factor of (k/4.5) 2 , where k is 
the polarizability of an atom and 4.5 is the numerical value of the polarizability of the hydrogen atom. 
In this way we can evaluate the binding energy for the state s = of the external electron of the ion of 
Cs~, at B — 30 T, as — Eq ~ 20 meV. For the same state of the ion of Ba~ we obtain, at the same field 
strength, — £o ~ 10 meV. As we see, the magnetically induced anion bound states for atoms with large 
polarizabilitics can very well be detectable in laboratories. 

Because of the large extension of the wave function of the loosely bound electron our results may also 
apply to molecules and clusters although in a less quantitative manner. As molecules and clusters can 
exhibit large polarizabilities |l7j we expect that magnetically induced bound states exist for them and 
possess considerable binding energies. 

One important issue which was not addressed in this paper are the finite nucleus mass effects which 
can affect the ion states significantly. The underlying physical picture of the finite nuclear motion for 
charged systems was discussed in detail in, e.g., Ipcfl . Here we only mention that the electron states are 
affected by an oscillating electric field introduced in the internal (electronic) system of the ion due to the 
rotation of the ion as a whole over a Landau- like orbit. For a crude estimate, we may assume that due to 
the ion motion the ion's states get the oscillator- like energy excess hflN, where f2 = eB/(Mc), M is the 
ion mass and N = 0, 1, 2, . . .. Since in the presence of the magnetic field the internal (electronic) degrees 
of freedom are inherently coupled to the center of mass ones [^(J , we can expect that similar to the case 
of the positive He + — ion |2l[] the minimal possible value of N for stable states coincides with s. In this 
case, an s— state of the external electron is still bound if — e s > hQs, or 

B 1.58 x 10 4 s 1 ' 2 , , 

Tt > ' (46) 

where A is the mass number of the ion. According to this criterion, for the H _ — ion the magnetic field 
strength which stabilizes the state s = 1 against the motional Stark effect, is large, B w 3.5 x 10 3 T, 
and cannot be achieved in laboratories. However, for heavier ions, the magnetic field which can stabilize 
them in excited states are smaller. For example, for the Cs~— ion, such a field is w Z.^^Sj 1 T, so the 
external electron can be bound in a few s— states at the magnetic fields available at laboratories. We must 
remark that these estimates of the center of mass effects on the ion bound states are very preliminary 
and based on the intuitive picture of the inclusion of the oscillator-like energies associated with the ion 
motion across the magnetic field in the total energy spectrum. We plan to perform the corresponding 
analysis accurately in order to see how the nuclear motion influences the induced bound states in the 
magnetic field. But already now we can speculate that the finite nuclear mass effects make at least the 
highly excited s— states unbound which means that in reality the actual number of bound ion states is 
finite and not infinite. 



11 



Acknowledgements 



We are grateful to M.V. Ivanov for providing us with numerically computed wave functions for the 
ground state of an electron moving in the Coulomb fields of two static charges. One of the authors 
(VGB) is pleased to acknowledge the Alexander-von-Humboldt Stiftung for the Fellowship. Support by 
the Deutsche Forchungemeinschaft is greatly acknowledged. 



Appendix 



According to Eq. ( |l0|) , the non-adiabatic correction to the potential of interaction between the hydrogen 
atom and an extra electron reads 

2 



hoo(r) 



2m r 



dri 



9 A I \ 



(Al) 



Here integration is performed over the position r\ of the core electron and the wave function <j>$ describes 
the ground state of the atom perturbed by the interaction with the static external electron. When the 
interaction is weak, this wave function can be evaluated by the first order perturbation theory, 



Mri,r) 



E 



W /0 (r) 



fj&O E 



(0) 



E 



(0) 



(A2) 



where 



and 



4^ are the wave functions of the ground and excited states of the non-perturbed atom, 



respectively, E^ and E^p are the corresponding eigenenergies and 

W f0 (r) = (<pf ) \W(r 1 ,r)\$ ) ) 
are the matrix elements of the perturbation operator 



e 
r 



(A3) 



(A4) 



Substituting (A2) into ( [AID and integrating over r% yields, due to the orthogonality of the hydrogenic 
wave functions, 



2m, 



E 



\ s fo\ 2 



(4 



.0)__g(0) 



(A5) 



where S — dW/dr. To simplify the latter expression, let us approximate the energy difference in the 
denominator by the quantity — Ry. Then, using the completeness property of the hydrogenic wave 
functions, we obtain 



/loo = 



2m c Ry 



(5 2 ) 00 - (SooY 



(A6) 



where (S 2 ) 0Q and Sqo are the mean values of the operators S 2 and 5, respectively, for the ground state of 
the non-perturbed hydrogen atom. To calculate them, it is most convenient to specify the internal atomic 
coordinates with the quantization axis along r. In these coordinates, assuming r > ri, we can replace 
the operator (A4) by the first, dipole, term of its multipole expansion (cf. Eq (|l5|)), W = (e 2 /r 2 )zi. 
Then we have S = — (2e 2 /r 3 )zi. For the latter operator, because of the z— parity of the hydrogenic state, 
one has Sqo = and Eq. (A.6) reduces to 
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2h 2 e 4 (zl) m 

hoo = -Vf^ • (A7) 

m c Ky r° 

Substituting there the known value (^i) 00 = «o an< ^ introducing atomic units, we obtain the estimate of 
the non-adiabatic correction given by Eq. (f44|), hoo{T) — 8/r 6 . 
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Figure Captions 

Figure 1. The one-particle potential of interaction between the static charge — e and the hydrogen 
atom (solid line) and the non-adiabatic correction to this potential for the cases when the charge is 
associated with the electron (long-dashed line) and with the muon (shot-dashed line). Dots show the 
reference data for the potential obtained from The part of the solid line connecting the dots is 

obtained by the spline interpolation, at lower r the solid line corresponds to Eq. (|4^) while at large r 
it corresponds to Eq. (|4l|) . Open circles show the non-adiabatic corrections obtained numerically and 
dashed lines correspond to the perturbation estimates. 

Figure 2. Binding energies for the states s = 0,1,2,3,4 (in meV) as functions of the magnetic field 
strength (in Tesla). The conversion rate used is: 1 a.u. = 2.35 x 10 5 T. Solid lines show the estimates 
for the ion H~ given by Eqs. (^) and ((3^), and dots represent the results of our numerical treatment. 
Dashed lines show the corresponding energies for the muonic ion, H/i~ . 
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